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Abstract 

We examine the time evolution of a quantized field in external backgrounds that violate the 
stability of vacuum (particle-creating backgrounds). Our purpose is to study the exact form of the 
final quantum state (the density operator at the final instant of time) that has emerged from a 
given arbitrary initial state (from a given arbitrary density operator at the initial time instant) in the 
course of evolution. We find a generating functional that allows one to obtain density operators for an 
arbitrary initial state. Averaging over states of the subsystem of antiparticles (particles), we obtain 
explicit forms of reduced density operators for the subsystem of particles (antiparticles). Analyzing 
one-particle correlation functions, we establish a one-to-one correspondence between these functions 
and the reduced density operators. It is shown that in the general case a presence of bosons (e.g., 
gluons) in the initial state increases the creation rate of the same type of bosons. We discuss the 
question (and its relation to the initial stage of quark-gluon plasma formation) whether a thermal 
form of one-particle distribution can appear even if the final state of the complete system is not in 
thermal equilibrium. In this respect, we discuss some cases when pair-creation by an electric-like 
field can mimic the one-particle thermal distribution. We apply our technics to some QFT problems 
in slowly varying electric-like backgrounds: electric, SU(3) chromoelectric, and metric. In particular, 
we analyze the time and temperature behavior of the mean numbers of created particles, provided 
that the effects of switching the external field on and off are negligible. It is demonstrated that 
at high temperatures and in slowly varying electric fields the rate of particle-creation is essentially 
time-dependent . 

Keywords: external particle-creating field; evolution of an arbitrary initial state; creation rate; 
electric, SU(3) chromoelectric and metric fields; initial thermal distribution. 



1 Introduction 

The effect of particle creation from vacuum by an external background (vacuum instability in external 
fields) ranks among the most intriguing nonlinear phenomena in quantum theory. Its theoretical analysis 
must be nonperturbative, and its observation in experiment is to verify the applicability of a theory 
in the domain of superstrong fields. The study of this effect started in connection with the so-called 
Klein [T] paradox, and was carried on by Schwinger [2] , who calculated the vacuum-to- vacuum transition 
probability in a constant electric field. A complete study of particle creation from vacuum by a constant 
electric field is presented in [31 S] . It should be noted that the effect can actually be observed as soon 
as the external field strength approaches the characteristic value (critical field) Ec — m?c'^ /\e\h ~ 1,3 • 
10"^^ V/cm. Although an actual possibility of creating these fields under laboratory conditions does 
not exist at present, e^e~-pair production by a slowly varying (external) electric field from vacuum is 
possibly relevant to phenomenology with the advent of a new laser technology, which may access the truly 
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strong-field domain. Electron-positron pairs can also be created by perfectly perturbative processes in 
crossed laser beams, when there occurs the merging of several photons. The control over the dependence 
between the power and frequency of lasers makes in possible to determine which of these mechanisms is 
responsible for pair-creation. This topic is widely discussed [5] at SLAC and TESLA X-ray laser facilities. 
Such strong fields may be essential in astrophysics, where characteristic values of the electromagnetic and 
gravitational fields in the vicinity of black holes are enormous. Electric fields near a cosmic string can 
become extremely strong [B]. In this respect, one needs to mention that the Coulomb field of superheavy 
nuclei may create electron-positron pairs; see [7]. Apart from purely QED problems, there arise problems 
of QFT in which vacuum instability in various external backgrounds plays an important role, for example, 
phase transitions in non-Abelian theories, the problem of boundary conditions, or the influence of topology 
on vacuum, the problem of a consistent vacuum construction in QCD and GUT, multiple particle-creation 
within the context of heavy- ion collisions, and so on. Particle creation by background metrics is important 
in black- hole physics [HI [S] and also in the study of the dynamics of the early Universe [T^ . Recently, it 
has also been recognized that the presence of a background electric field must be taken into account in 
string theory constructions; see, e.g., [11: and references therein. 

Considerable attention has been recently focused on a non-perturbative parton production from vac- 
uum by a classical chromoelectric field of SU(3) [I2l[l3] and SU(2) [M], in the framework of a modern 
version of the known chromoelectric flux tube model (TS], the latter being an effective model for the 
confinement of quarks in QCD (previously pair-creation by a constant fleld has been calculated for SU(2) 
in [TO] and for SU(3) in XU). The model ensures a very good description of the phenomenology of hadron 
jets in high-energy — e~ and p — p collision experiments (a further development of the basic model 
and phenomenological applications can be found in the review 18J). This model probably describes the 
initial stage of quark-giuon plasma formation reasonably well (in particular, the transversal spectrum 
of produced soft partons). Such a state may be produced at high-energy large-hadron colliders such as 
RHIC (Au-Au colhsions at Vs = 200 GeV) [19] and LHC (Pb-Pb colhsions at Vs = 5, 5 TeV) At 
present, this initial stage is related to an effective theory, the color glass condensate [2l] (see also the 
review papers |22j). which is the coherent limit of quark-gluon plasma at high energies. In such a picture, 
after a nuclei collision, a strong classical chromo-electric-magnetic field is created due to relatively slow 
fluctuations of color density. Such a field is sufficiently uniform in the direction that is transversal to 
the beam direction and has a longitudinal chromoelectric component |14[ 123] . This component is much 
more intensive than the transversal component; see [14 . Thus, the color glass condensate picture pro- 
vides strong arguments in favour of the chromoelectric flux tube model and allows one to calculate field 
configurations in a tube. In particular, such a physical picture allows one to accept a quasi-constant 
chromoelectric field as a good approximation at the above-mentioned initial stage. It should be noted 
that experimental data on heavy-ion collisions that exist at present can be interpreted as quantum parton 
production by an external chromoelectric field from vacuum and many-particle states. 

There is a considerable interest in particle-creation at finite temperatures and at a finite particle 
density, which is basically motivated by heavy-ion collisions, cosmological QCD phase transitions, and 
dark matter formation. For example, thermally-influenced pair-production in constant electric flelds at 
the one-loop level has been searched for in [24l [25l [26l [27] . 

The above calculations have been carried out within the theory of a quantized field placed in an 
external background. A consistent description of a complete QED (interacting quantum electromagnetic 
and matter fields in a particle-creating background) with unstable vacuum that treats the interaction 
with external backgrounds nonperturbatively has been developed in [28 . Possible generalizations of the 
formalism to external gravitational and non-Abelian gauge fields have been presented in [5S1 [30] ; see also 
[3T] . An attempt to extend this technics to the thermal case has been taken in [32] . Calculating particle 
creation by black-hole metric. Hawking has discovered that the density matrix of created particles at 
spatial infinity has a thermal character. The question arises, is such a character related to the particle- 
creation mechanism in general, or to the gravitational origin of the background? A way to answer this 
question is to elaborate an adequate technique which could allow one to include arbitrary mixed initial 
states, in particular, thermal initial states, in the corresponding particle creation formalism [28] . 

In this article, we present a development of the particle-creation formalism [28] that is capable to 
answer some of the above questions. We examine the time evolution of a quantized field (bosonic or 
fermionic) in external backgrounds that violates the stability of vacuum. In fact, we deal with a quadratic 
field theory of noninteracting (between themselves) particles. Our purpose is to study the exact form of 
the final quantum state (the density operator at the final instant of time) that has emerged from a given 
arbitrary initial state (from a given arbitrary density operator at the initial time instant) in the course 
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of evolution. 

The article is organized as follows. Section 2 has an original, albeit rater technical character. There, 
we derive exact expressions for density operators (more appropriate for the generating density operator) 
by applying the path integration method. Some necessary formulas are placed in Appendix. In Section 
3, having an exact expression for the generating density operator, we obtain reduced density operators 
for the subsystems of particles and antiparticles. We introduce and calculate the one-particle correlation 
functions and establish a one-to-one correspondence between these functions and the reduced density 
operators. In particular, this allows us to restore the reduced density operator of the complete system 
starting from the one-particle distributions (of course, this is possible only in the model under consider- 
ation being a quadratic theory). It is demonstrated that in the general case the presence of bosons (e.g., 
gluons) in the initial state increases the creation of the same kind of bosons. We discuss in detail the 
question (and its relation to the initial stage of quark-gluon plasma formation) whether the thermal form 
of the one-particle distribution can appear even if the final state of the complete system is not in thermal 
equilibrium. In Section 4, we discuss the expressions we obtain for the density operators and one-particle 
distributions in electric-like backgrounds: electric, SU(3) chromoelectric, and metric. In particular, we 
analyze the density operators and one-particle distributions in the so-called T-constant electric back- 
ground (a field exists during a finite period of time T), and demonstrate how such a problem is related to 
particle creation by gravitational fields (Hawking's effect). We present some examples when pair-creation 
by an electric-like field can mimic the one-particle thermal distribution. Then, we analyze the time and 
temperature behavior of particle-creation when the effects of switching on and off are negligible. In par- 
ticular, we demonstrate that at high temperatures the production rate is non-trivially time-dependent. 
This result has to be taken into account at high temperatures. 

2 Density operator in pair-creating backgrounds 

We consider a quantum field ^'{x) in an external background. The quantum field can be scalar, spinor, 
etc., and the background can be the classical electromagnetic, Yang-Mills, or gravitational field. In the 
general case, the background is intense, time-dependent, and violates the stability of vacuum. Such a 
background must be treated nonperturbatively. We plan to follow the formalism proposed in !28| . 

2.1 Some relevant relations 

It is assumed that there exists a set of creation and annihilation operators ajj(ti„), a„(ii„) of particles 
o-n{tin), 0'n{tin)i ^nd antiparticlcs 6jj(ii„), bn{tin), respectively, at the initial time instant ii„ (ij„ — > — cxd), 
and a set of creation and annihilation operators of particles, a]i(ioMt)i CLnitout), and antiparticles, &Ji(iojii)i 
bn{tout), at the final time instant tout {tout — *■ oo). By n we denote a complete set of possible quantum 
numbers. The total Hamiltonian H (t) of the quantized field under consideration is diagonalized (and 
has a canonical form) in terms of the first set at the initial time instant, and is diagonalized (and has 
a canonical form) in terms of the second set at the final time instant. Nonzero (anti)commutator|3 are 
given by 

[o-n (^m ± — [an(io«t)' '^m(^ojii)]± 

= [bn{tin),b\^{tin)]± = [6„ (tout ) ' ^rn (*o«t )] ± = ^nm ■ (1) 

The vacuum states |0,ti„) at tin and \Q,tout) at tout are defined as usual: 

a{tin)\Q,tin) = 6(t„)|0,t„j) = 0, a{tout)\^,tout) = b{tout)\^,tout) = 0. 

Therefore, we define these states as the states that minimize the mean value of the Hamiltonian H (t) 
at t = tin and t = tout, respectively. According to the general princilples of quantum theory, the 
probability amplitude of transiltion from a certain initial state bl^{tin) ■ ■ ■ ajj(ii„)|0, tin) to a certain final 
state bl^,{tout) ■ ■ ■ o-li' (tout) \0, tout) in the Schrodinger picture has the form 

{O, tout I O^n' {tout ) • ■ ■ bni' {tout)U {tout, tin) bln{tm) ■ ■ ■ aji(tm)| Oi ^m) , 

^The subscript — denotes the commutator (in the Bose case), whereas the subscript -|- denotes the anticommutator (in 
the Fermi case) 
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where U {t, t') is a unitary evolution operator of the system. Let p {tin) = P (aU^m), o-{tin), b'^ {tin), b{tin)) , 
tip {tin) = 1, be the density operator of the system under consideration at the initial time instant. As 
the system develops in time, this density operator becomes p{tout) for antiparticles at the final instant 
of time: 

P {tout) — U {tout, tin) P {tin) U'^ {tout, tin) ■ (2) 

To pass to the Heisenberg picture, we introduce finite-time evolution operators f^(±), 

= U (0, t„0 , = U (0, tout) , U {tout,t„,) - ■ 

We then define a set of creation and annihilation operators a]j(m), an{in) of m-particles, as well as 
similar operators bl^{in), bn{in) of in-antiparticles, a corresponding m-vacuum |0, in), a set of creation 
and annihilation operators ajj, a„, of out-particles and similar operators bl^, bn of owt-antiparticles, and 
a corresponding out- vacuum |0), 

{a\in),a{in),b'^{in),b{in)) = ri(+) {a) {tin).,a{tin).,b\tin).,b{tin)) f^j^^ , 
{a\ a, b\ b) = {a'' {tout), a{tout), b'' {tout), b {tout)) , 

|0, in) = r!(+) |0, t^n) , |0) = r!(_) |0, tout) . (3) 

The in- and oui-operators obey the canonical commutation relations ([T]). 

The entire information concerning the processes of particle creation, annihilation and scattering is 
contained in the elementary probability amplitudes 

^(+l+)mn = c^^(0|a™at(m)|0,m), 
(-|-)mn = c;;^(0 \bnibi{in)\ 0,in) , 
«^ (0| - +)nrn = ^^7^(0 | &,\ (m)a]„ (m) | 0,m) , 

«^(+-|0)™„-c;i(0|a™6„|0,*7^), (4) 

where Cy is the vacuum-to-vacuum transition amplitude, Cy — {0\0, in) . 

The sets of in and oui-operators are related to each other by a linear canonical transformation (it is 
sometimes called the Bogolyubov transformation). As has been demonstrated in the general case, such a 
relation has the form (see [28] ) 

V{a\a,b\b) ^ ( 

a'^ {in) , a{in) , b^ {in) , b{in)) , 
\0,in) = V\0) {cu^{0\V\0)) , (5) 

where a unitary operator V has the form 

V = V4V3V2V1 , (6) 

an41 

Vi = expl—bw (0| — \-) a} , V2 = exp lnu> (+|+) a} , 

v:i — exp {—nblnw (— |— ) b^} , V4 ~ exp { — na^w (H — |0) 6^} , 



1 Fermi case 
— 1 Bose case 



(7) 



Using this expression for V, one can find 

Cy = {0\V\0) =exp{-Ktrlnu;(-|-)} . (8) 
The density operator p of the system under consideration in the Heisenberg picture is defined as 

p = il(^+)p{t^n)^l_^^) ^ p{a^{in),a{in),b^{in),b{in)) , 

p = n(^_)p{tout)^]),irp^l. (9) 



^We use condenced notation, for example, 

bw (0| h) a = ^ b„w (0| \-)nm ■ 
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Suppose a physical quantity is given by an operator F {tout) at the final time instant as 

F {tout) = F {a^tout),a{tout), b\t out), b {tout)) ■ (10) 
Then its mean value at the final instant of time is given by 

(F) = tr \f {tout) P {tout)] = tr [ Fp] , (11) 

where 

F = r!(_)F {tout) ^\_) = F (at, a, b\ b) (12) 

is the operator of a physical quantity in the Heisenberg representatiorH. 

The amplitudes (|4]) can be calculated with the help of some appropriate sets of solutions of the 
corresponding relativistic wave equation (RWE) with an external field (Klein-Gordon, Dirac, linearized 
Yang-Mills), see [2Hll3Q], as follows. Namely, the m-particles are associated with a complete set (in- 
set) of solutions of the RWE {c^^'ni^)} with asymptotics i^ipj^{tin,^) at the initial time-instant tin being 
eigenvectors of the corresponding one-particle Hamiltonian Ti.{t), 

'H{t,n)iiJn{Un,x) = C£i^\i^n{ttn,^) , (13) 

where Sn"^ are the energies of in-particles in a state specified by a complete set of quantum numbers 
n, and el^^ > 0. Here and elsewhere, ^ = ±, being (-I-) for particles and (— ) for antiparticlcs. For 
the sake of simplicity of exposition, we use such a gauge-fixing of the external field that condition (|13p 
should have the same form as it has in the absence of the external field. The owi-particles are associated 
with the complete out-set of solutions {^■(/'n {^)} of the RWE with asymptotics ''^'^{tout,^) at tout being 
eigenvectors of the Hamiltonian Ti.{t) at tout, namely, 

n{tout) ^^n{tout, X) = al^^ ^V„(iont, X) , (14) 

where e^^-* are the energies of owi-particles in a state specified by a complete set of quantum numbers n, 
and el^^ > 0. We suppose that the external field is such that soluitons of this eigevalue problem actually 
exist. 

The out-set can be decomposed in the in-set as follows: 

ii,{x) = +^{x)G (+1^) + K^^{x)G , (15) 

where k = +1 for fermions and k = — 1 for bosons. The decomposition coefficients G ^(jj'- ^ are expressed 
via inner products of these sets, G ( |^ j — ( cV'mj^ ''Pn ) • These coefficients obey the unitarity relations 

G (^|+)G(+|';)+«G(«|_)G ^C"^, 

G(cl+)G(+|c)+«G(cr)G(-|c)=c'^, 
G(+|+)G(+|-)+«:G(+r)G(-|-)=0, 

G{+\+)G{+\-)+nG{+\^)G{^\-)^0, (16) 



(<i<)^«(d<r 



that follow from the normalization conditions for the solutions. Here, the notation G | 
has been used. 

The quantum Heisenberg field 4'{x) can be expressed both from the creation and annihilation operators 
of in-particles and from the creation and annihilation operators of oui-particles: 

'ij'{x) = -\.ij}{x)a{in) + ^tp{x)b'' {in) = '^ijj{x)a + ^ip{x)b^ . 

This makes it possible to express the linear canonical transformation between the sets of in and out- 
operators in terms of the coefficients G (^^ , 

a = G{+\+)a{in) + G{+\-)b\in), 
nb^ = G{-\+)a{in) + G{-\-)b^in), (17) 



^All operators in the Heisenberg representation are denoted by the turned over hat in what follows, e.g. A. 
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and their Hermitian conjugated ones. Using relations PT|) . one finds the amphtudes Q in the form: 



w 


(+I+) = G(+|+ 








w 


(0|-+) = -G 


'-\-y'G{-\+)=KG 






w 


(+ - |0) = kG ( 


+\+y'G{+\-) = -G 






the same relations, 


one finds that 


the differential mean 


numbers 





(18) 

, one finds that the differential mean numbers Hm'' of particles/antiparticle 
created from vacuum are 

= (0,zn|aUn|0,m)= [G(+|-)G(-|+)]„^^ , 

= {0,zn\blb^\0,tn)=[G{-\yG{+\-)]^^. (19) 

Using the unitarity relations (|16p and representation ([8]), we find that the probability of a vacuum to 
remain a vacuum is given by 

= Icp = exp{Ktrln[l-AiG(-|+)G(+|-)]} . (20) 

Therefore, we can see that the non-stability of vacuum in external fields is mainfest in the fact that 
P" < 1, which is in one-to-one-correspondence with the fact that the differential mean numbers are 
different from zero. We can now establish the following general property. Let the time-dependent external 
field turn on at the time instant ti and turn off at the time instant and in general behaves arbitrarily, 
on condition that the solutions of problems and ((M]) do exist. Then, the in-set of solutions {c'0ti(^)} 
is defined equally for any solution of the time instant tin < ti, and the out-set of solutions {^V'n i^)} 

is defined equally for any choice of the time instant tout > ^2- Therefore, the coefficients G |<^j, and 
consequently all the above-calculated amplitudes, depend only on the dependence of the external field at 
the finite interval of time between ti and t2- This also implies that the in- and out- vacua, as well as the 
in- and out-operators of creation and annihilation, do not depend on the choice of the time instant tin 
and tout, respectively, if tin < ti and tout > ^2. 

2.2 Generating density operator 

We introduce the following generating operator R{J): 

R{J)^^R{J) , trP(J) = l, (21) 
R{J) =M„exp a) {in) - l) a{in) + b\in) (j^"' - l) h{in) 

where the Grassmann-even variables J = (^•^n ^^ ^'^''^ sources; JJmn = SmnJn ''', A/in is the sign of the 

normal form with respect to the in- vacuum, and Z — trP(J) is a normalization factor (statistical sum). 

In order to complete the calculation, it is efficient to use the path integral representation. In the 
fermion case, we use a path integral over anticommuting (Grassmann) variables, which is understood as 
Berezin's integral [33] at k — 1, 



. ^ — Ka^Ka 



:= det if" : J exp {k\*K-^\ + aU + \*a} Ud\*d\ :, (22) 



where a^, a are some creation and annihilation operators, and : . . . : realizes the normal form of the 
operators a\ a. All operators and a can be considered as Grassmann-odd variables under the normal 
form; therefore, we can calculate the complete path integral as a Gaussian integral over Grassmann- 
odd variables. In the boson case, we use the path integral over commuting variables at k = — 1. In 
this case, we can regard all the operators and a as bosonic (ordinary) variables under the normal form 
sign, so that the path integral in ([22]) is a usual Gaussian path integral, where X*K~^X > 0. 

Note that the trace of the normal product of creation and annihilation operators can be calculated 
by using the path integral representation, according to ()110|) . see Appendix. For example, by calculating 
Z we obtain 

Z = expi AtV [infl + KjjWjj [ln(l + Kj(-))j \. (23) 

1 ^ — ^ L \ /J nn ^ — ^ L \ / J mm 
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Having at our disposal the generating operator (|2ip . we can obtain the different density operators (in 
the Heisenberg representation) corresponding to different initial states of the system. We represent some 
examples below: 

a) By setting J = 0, we obtain a density operator of the system that is found in the pure vacuum 
state at the initial time instant, 

p„=i?(0) . 

Indeed, using relation (jlOSp from Appendix, we obtain 



=A/i„exp|— {in)a{in) + {in)h{in)\^ — |0,m)(0,m| 
In addition, we define the following generating functional of moments: 



(24) 



(J) = (0, in\ exp 



|0,m) = iY ((>{,]) 



4> (J) = exp 



Pv 



(25) 



which is useful to investigate the final state evolved from vacuum at the initial time instant. 

b) The density operator P{rn}M-{n}N '-'^ system which is found in a pure state with M particles 
and N antiparticles (with the quantum numbers {mi, . . . ,m]\i} = {m}]\i and {ni, . . . , rt^v} — {j^lw, 
respectively) at the initial time instant can be obtained from the generating operator R{J) as follows: 



P{m}M;{n}f 



= l*{m}M;{n}„ («'^)>(*{m}M;{n}„ (*") I , 



(26) 



j=o 



where 



M AT 

' ;{n}iv 

1=1 3 = 1 

N M 

{in) I = (0, ot| ]^ bn. {in) J]^ am, (in) 

3 = 1 i=l 



c) Let us set 



(27) 



where ei^'' are the energies of particles or antiparticles with the quantum numbers n; /i^^-' are the corre- 
sponding chemical potentials, and O is the absolute temperature. One can see that with such a choice 
of sources the generating operator (PT|) becomes the density operator of the system that has been in 
thermal equilibrium at the initial time instant. Using relation (jlOSp from Appendix, we obtain an explicit 
expression for p^, namely. 



Pp^ R\e 
Z — exp 



= ^cxp 



|- a){in)E^+^a{in)+b\in)E^-^h{in) | 



L^ln(l 

I n 



+ H 



1 + KC 



or 



Pfi^ Z ^ exp ^ -/3 



(28) 



where H is the Hamiltonian of the system (written in terms of in-operators); iV^^^ are the operators of 
m-particle or m-antiparticle numbers, 

H = a^in)e^^^a{in) + b^in)e'^^^ b{in) , 
iV(+) = at(m)a(m) , = b\in)b{in) , 

and the matrices -E^^^ and e'''^ are defines as Emn = SmnEn ^ and £ml = Smn^n^- 



7 



We can see that the problem of calculating the mean value of an operator F {tout) for the system 
being at the final time instant is related to the problem of calculating the quantity tr [FR (J)] , where 
F is a. Heisenberg operator corresponding to F {tout)- Such a quantity can be represented as follows: 

oo 

tr[FR{J)]^ E E J^i'^i{^}MAn}j,)\FR{J)\^{{m},,,{n}j,)), 

M,N=0{m}{n} 

i{m}M , {n}N)\ = (0| bnr, ■ ■ ■ hn^a^n^ ■ ■ ■ a,m (29) 

Calculating tr [FR{J)'\ according to (|29p . it is convenient to have an expression for the operator R{J) 
in terms of out-operators. One can see that such an expression has the form 

R{J)^VU{J)V\ C/ ( J) =: exp at - l) a + 6t (^j(-) - l) 6 :, (30) 

where : • • • : is the sign of the normal form with respect to the out-vacuum, and the operator V is defined 
by (O. A normal form of the operator R{J) with respect to the out- vacuum is calculated below. 

2.3 The normal form of the generating operator 

First of all, we rewrite the operator expression ([30|) as follows: 



R{J) = ^iV3V2Y{J)vlvlvl, Y = v,U{J)vl, (31) 

where the operators Vi, i = 1, ...,4, are given by ([7]). Using formula l|104p from Appendix, we represent 
the operator Y{J) in the form 

Y{J) = Y (J) U{J) ,Y (J) = exp {~bBa) exp (-aU (J) b'') , 
A{J)=I^+^B^^-\ B = w{0\-+) . (32) 

Both operator exponents in the expression for Y (J) can be written in terms of Gaussian path integrals. 

Consider, first of all, the fermi-particle case. In this case, we can treat the anticommuting operators 
a and b (or a+ and 6+) as Grassmann-odd variables. Then, according to representation (|22p at k = 1, 
we have 

y = det A det B / exp f A*B"^ A + A* a) <^Ud~X* d~XdX* dX , 



$ = exp [bX + X aj exp (aU + X*b^) . (33) 

With the help of relation p06p from Appendix, we represent the operator $ in the normal form 

$ =: exp (aU -I- A*6t -|- 6A -I- X*a + X* X + AA*) : . 

Then, using formula we can calculate the complete path integral In the Bose case, we can 

examine all the operators a\ b\ a, and b as bosonic (ordinary) variables under the normal form sign, so 
that the operator Y (J) can be represented as a usual Gaussian path integral by applying representation 
(|22)) at K = —1. Calculating these Gaussian integrals, we obtain the normal form of the operator Y, 

Y = det (1 + kAB)"^ : exp 

A++ = kAB (1 + kAB)^^ , A^_ = kBA (1 + kBA)^^ , 

A+- = (1 + kAB)^^ a , A-+ ^ B{1 + kAB)^^ . (34) 
Using relation (jl07p from Appendix, we represent the operator Y in the normal form 

Y = dot (1 + kAB)'' : exp |-aU++a - b^ A—b - atl+_&t „ bA_+aj : , 
I++ = 1- (1- A++)J(+), = 1-(1-A__)J(-), 

I+_ = I-+ = J^"U-+J(+^ . (35) 
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With the help of relation ([5]), we rewrite the operator V3 as follows: 

U3 = exp [— K61nu; (— |— ) 6^] = c„ exp 6^ In w (— | — 6 

Then, using formulas (|105p derived in Appendix, we represent the operators V3V2 and vlvl from pTjl in 
the normal form as follows: 



V3V2 = Cy : exp 

t t * 
"2^^3 = : exp 



(w {~\^f - 1) &] exp [flt (+|+) - 1) a] : , 
(+1+)^ - 1) a] exp [fo^ (w H- f^ - l) & 

Finally, using relation (|107p . we obtain the normal form of the operator R (J) , 

R{J) = |c„|^det(l + KAB)'' :exp [-a'' {1 - D+) a - b^l - D^) b - C%'' - bCa] :, 
D+^w (+1+) (1 + KABy^ JJ^+'w (+1+)^ , 
= u;(-|-)tj(-)(i + KBA)"'u;(-|-) , 
C^w (-l-)t (1 + nABy^ j(+)u> (+|+)t + (+ - |0)^ . (36) 

Representation is useful, since it allows one to calculate the trace (I^H]) by using the path integral 
techniques described in Appendix; see, eq. (|110p . 

As an example, let us consider the density operator defined by ([24|l . Using ([36|l . we represent this 
operator in terms of owi-operators, as well as in the normal form 



p,^ = R{0)^ |c„|^ : exp [-a^a- 5^6- Ka^w(+- |0) 6^ - nbw {+ - |0)'^ a 



(37) 



In a similar way, the operator (j){J) in the expression of generating functional ()25p can be transformed 
to the normal form 



(J) = |c„|2 : exp [-a^a - b^b ~ Ka''e^''^^w {+ - |0) e^' ' b^ - Kbw {+ - |0)^ < 



(38) 



where the representation (|37|) for p„ has been used. Then, applying the path integral representation for 
traces (jllOp . and using formula (|107p . we represent the generating functional of momenta as follows: 



(J) = |c,|'exp|Ktrln [l + (+ - |0)t e^'^' w (+ - |0) e^'^'j } 



(39) 



3 Reduced density operators and correlation functions 

3.1 Reduced density operators 

In the general case, the states of the system under consideration at the final time instant contain both 
particles and antiparticles due to the pair-creation by external fields and the structure of the initial state. 
On the other hand, we are often interested in physical quantities F± which describe only particles (+) or 
antiparticles (— ) at the final time instant. The corresponding operators F± are functions of either a^,a 
or 6^,6, 

F+=F+{a\a), F^^F^{b\b). (40) 

The mean values of operators F± and the entire information concerning the subsystems of particles and 
antiparticles can be obtained from the so-called reduced density operators, which we shall define below. 
We present the basis vectors from (1^^ as follows: 

I* ({m}^, , {n}^)) - 1*, ({m},,,)) (g> \^f, ({n}^)), |0) = |0), (g> \0), , 

l*a ({m},,)> - al^ . . . al^,^ |0), , {{n}^)) = bl . . . fot „ |0), , (41) 

where |0)a and |0)b are the vacuum vectors of particle and antiparticle subsystems. The mean values of 
the operators F± are 

(^^±) - tr+tr. (F±p) , (42) 
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where ~p is the density operator of the system, and the reduced traces tr± of an operator A are defined as 



tv+A 



M=0 {m} 



*a ({m}^,) lAlM/, ({m}^)) 



tr_i = ^ ^ (Af ({mU,) ({m}^^)) . (43) 

A/=0 {m} 

We define the reduced density operators (in the Heisenberg picture) of the subsystems of particles 
and antiparticles, respectively, as follows: 

P±=tT^p. (44) 
Then mean values (|^^ can be calculated with the help of the reduced density operators p^ as follows: 



(^^±) = tr± (F±p± ; 



(45) 



Even if the initial state of the system is a pure state, the reduced density operators pj. describe mixed 
states. In some physical problems, the use of a reduced density operators is inevitable. For example, 
considering particle creation by the gravitation field of a black hole we have only the reduced operator 
of the particles created outside the black hole, since we do not have any information about the particles 
behind the horizon, [S] [9] . 

In a similar manner, we introduce the reduced generating operators R± (J) as follows: 

R± (J) = tr^R{J) . 

Using the path integral representation for traces (|110p . representation as well as p07p . we obtairQ 

R+ (J) = Z^^ : expj-a^ (1 - K+ (J)) a} :, 
R- (J) = ZZ^ : cxp {-6t (1 - (J)) b} : , 



K± (J) = D±+C^ {1 + " C , 

(J) = Z-'^ \cX det (1 + nABf det (1 + nD^f 



(46) 



The reduced generating operators R± (J) allow one to obtain the reduced density operators pj. for 
different initial states of the system. Consider below some examples: 

a) Selecting all J = in (|46p . we obtain reduced density operators p^^ = R± (0) of a system that has 
been in a pure vacuum state at the initial time instant. Explicit expressions for R± (0) follow from (|46p 
with account taken of 

K± (0) ^wi+-\0)w{+~ |0)t , (0) = |c„|' . (47) 

The same result has been obtained in [531 131] by a straightforward calculation. 

b) The reduced density operators Po-n± and Pm-o± of ^ system that has been in a pure state with 
particles or antiparticles, respectively, at the initial time instant can be obtained from the generating 
operator_R^ (J) = ZR± (J) as follows: 



Pm:0- 



dR+ (J) 



(+) 



9J, 
dR_ (J) 



\a^w ( 



{-) 



PO;m + 



dJ, 
dR+ {J) 



,7=0 



J=0 



[^i-\-)b']„/P.- 



bw{-\-y 



dJ, 



(-) 



Pv^ 



w{-\-)w{-\-f 



J=0 

a^w(+- |0)i«(-|-)t 



Pv+ 



w{-\-)w{+~~ |0)^a 



Pm:Q - 



dR_ (J) 



dJ, 



(+) 



= Pv- 



W {+\ + )^ W {+\+) 



.7=0 



[b^w (+ - |0) w (+1+)*] ^ p,_ w {+\+f w{+~ |0)^ b 



*It should be noted that the symbols of the normal form of the operator have been presented in 1241 via some path 



integrals. The explicit form of operator has been presented for J, 



(C) 



Unfortunately, it contains some misprints. 
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c) Let us set the sources in as in (|27p . One can see that for such a choice of sources the reduced 
generating operators (|2ip become the reduced density operators of the system that has been in 
thermal equihbrium at the initial time instant. 



3.2 One-particle correlation functions 

Let us now examine the following generating functions: 

= tr {alia„iR) = tr+ (al^a„iR+) , 
nU = tr {bib„,R) = tr_ {bib„,R+) . (48) 

They generate one-particle correlation functions for different initial states of the system. Setting the 
sources (taking the corresponding derivatives, if necessary) in (|48p as has been demonstrated in Sec. 2, 
we choose the required initial states. The diagonal elements NmL are the generating functionals for the 
mean numbers Nm"^ of particles/antiparticles with quantum numbers m at the final time instant (further 
differential mean numbers). In what follows, we refer to quantities (I48p as correlation functions. 

The correlation functions Nnm can be expressed via the matrices K^^ and vice- versa, as follows: 

Note that the quantities are functions of elementary probability amplitudes (|4|) . 

Relations (|49p can be proved as follows: first, using the commutation relations H]) we represent (|48| 
as traces of operators in the normal form: 

Ni+2 = tr+ [aiR+ {K+a}J , N^-J = tr_ [&t i?_ {K^b)J . (50) 
The quantities Nnm can be obtained from the generating functions (j, j) as follows: 



djndjr. 

where 



(51) 



Z+ (j, j) = Z-hr+ : exp {-at [1 _ IK+] a] : , 
Z_ (j,j) = ^lHr_ :cxp{-6t[l-Ii^_]fe} 

(0, 0) 1 , Iran = S,nn + Jmjn , C = ± , (52) 

and j and j are some new sources. Traces in (j52l) can be calculated by formula (jllOp from Appendix. 
Thus, we obtain 



Zc (j, j) = exp I « [In (1 + «Kc)]„„ I 



_ (53) 

Then, relations ([49]) follow from (l5l|) and ((53|) . 

The normalization conditions and the second relation (|49p imply that the quantities Zi^ can be ex- 
pressed in terms of N^''-' as 

Zc-exp|K^[ln(l + Acifc)]„„| =exp|-A^^[ln(l-.jN(^)^)]^J . (54) 

Now, we are going to relate the quantities Nnm with the correlation functions Nnm (in) of in-operators, 
Ni+2{tn) = tr [ai{tn)a^{tn)R] , NiJ{tn) = tr [bUtn)b^{in)R] . (55) 
Using representation (pij) for ^, one can see that 

tK) 

N(£(w) = Sr,rnN,^p [m), TVi^) (*,i) = ^ , (56) 

1 -\- hiJm 
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where Nm\in) are the differential mean numbers (generating functions for differential mean numbers). 
Indeed, let us take expressions (|48)) for N*^""^ via traces in the complete Fock space. These traces can be 
written in the m-basis {{m}j^ , {n}jY ! *")) = 1^* , {n}j^)). Using the canonical transforma- 

tion (I17p . we express the operators a^, a, 6^,6 via the operators a^in) , a{in) , b'' (in) , b{in) and calculate 
the traces explicitly. Then we obtain 

N(+)^ = G (+1+) N^+Hm)G (+1+) + G (+|_) [l - ^^(-^(in)] G (_|+) , 

N(-) = G (-|_) E^-Htn)G (_|-) + G (-|+) [l - KN(+)(m)l G (+|-) . (57) 



Thus, due to relations (|49)) . ([57]), we have explicit expressions for the complete generating density 
operator (|36|) and reduced generating density operators (j46|) via both correlation functions of m-particles 
and oui-particles, and via elementary probability amplitudes (jl]) as well. 

We stress that there is a one-to-one correspondence between the one-particle correlation functions, 
and the form of the reduced density operator of the complete system. This correspondence is related to 
the choice of a model, which is a quantized field placed in an external background. In fact, we deal with a 
quadratic system of noninteracting (between themselves) particles. Of course, this fact is well-known for 
free-particle systems. Our consideration generalizes this desription to the presence of a particle-creating 
background. For systems of interacting particles, there remains an important question: suppose the 
one-particle distribution at the final time instant is a thermal one. Can one assert that the complete 
system is in a thermal state with a given temperature (one that determines a one-particle distribution)? 
Such a question seems to be relevant to the problem of particle creation by black-hole gravitational fields 
(Hawking's radiation), where one-particle distributions of created particles have a thermal form. 

Below, we examine some illustrations of the previously obtained general formulas. 

Let the initial state of the system be vacuum (J = 0), then ([57]) reproduces formulas (fT9|) for the 

differential mean numbers Hm-* = Nm^ of particles/antiparticles created from vacuum by an external 

,7=0 

field. 

Let us examine a common case (for example, a uniform external field) when particle/antiparticle states 
are specified by quantum numbers (the same being valid for particles and antiparticles) that are integrals 

of motion. In this case, all the matrices G ^<^|^ ^ in (fT5)) are diagonal and the differential mean numbers 

pj)) of particles/antiparticles created from vacuum coincide: = Hm' = H„i. Using formulas (|57|) . 

(|19p . and the unitarity relations (jl6[) . one can obtain the expressions for the differential mean numbers 
of particles/antiparticles, namely, 

iVK) = (1 - kH„) 7V«) (m) + H„, [l - KiV,(-?) (m)] . (58) 

If the initial state differs from vacuum, the differential mean numbers of particles/antiparticles created 
by the external field are given by the difference Aivii^ = N^^ - N^\in). One can see that 



A7Vi+) = ATVi") = AN, 



m 1 



1 



K[N'i+\in)+N'^\in))\. (59) 



Even if H„j ^ 0, no particle creation of fermions with quantum numbers m occurs iiN^\in)^NL \in) = 
1. Since k = — 1 for bosons, is always positive and is larger than H^. That is, the presence of 

matter at the initial state increases the mean number of created bosons. 

In some articles devoted to the chromoelectric flux tube model (see, e.g., [12l [HI ES] ) , one encounters 
an (inexact) interpretation of the well-known Schwinger formulas describing pair-creation from vacuum 
by a constant electric field [5]. This interpretation may lead to incorrect results for some field strengths, 
as noted in [36]. Below, we discuss this problem and present correct relations that will be used in 
the subsequent section. We recall that, by using the proper-time method, Schwinger calculated the 
one-loop effective Lagrangian L in electric field and assumed that the probability P"" of no actual pair- 
creation occurring in the history of the field during the time T in the volume V can be presented 
as P" = |c„|^ = exp{— Fr2ImL} (for a subsequent development, see the review [27] )■ Schwinger 
interpreted 2ImL as the probability, per time unit, and per volume unit, of creating a pair by a constant 
electric field. Some arguments in favour of such an interpretation can be found, for example, in the 
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book [35] and the article (TS]. The interpretation remains approximately valid as long as the WKB 
calculation is applicable, that is, VT2lmL ^ 1. Then the total probability of pair-creation reads as 
1 — P" w VT2 Im L . To calculate the differential probabilities of pair-creation with quantum numbers 
771 (for instance, momentum and spin polarization), one can represent the probability as an infinite 
product: 

P'' = ]^e-2i«^^", (60) 

m 

where a certain discretization scheme is used, so that the effective action S — VTL is written as 5 = 
X^m '^rn- AH this is possible only if m are selected as integrals of motion. Then, e"^^™'^™ is the vacuum- 
persistence probability in a cell of the space of quantum numbers m. Using the WKB approximation 
in the case 2Im5'm «C 1, one obtains for the probability Pm of a single pair-production with quantum 
numbers m and for the corresponding mean values of created pairs the following relation: 

Prn~2lmSm. (61) 

By analogy with one-particle quantum mechanics, one usually rewrites ()61|) for fermions, 

H„, w-ln(l-P„0 w2Im5,„. (62) 

It is clear that and ([5^ coincide in the first order with respect to P„i . Then, it follows from ([50)1 
that 

P'^-Uil Pm) . (63) 

m 

Using the same analogy for bosons and rewriting (|6ip as 

H„«ln(l + P„)«2Im5™, (64) 
one obtains the following approximate relation: 

P'' ^l[{l+P,ny' . (65) 



It turns out that for the field under consideration, by using the WKB calculations and relations 
(|65| . one can reproduce Schwinger's result for P". This fact brings the temptation to interpret the latter 
formulas as exact ones, replacing there "~" by . However, one should say that such an interpretation 
is, in particular, equivalent to the assumption H.,„ = 2Im5'm. Nevertheless, as we shall demonstrate 
below, the latter relation is not an exact one and it is valid only in the approximation 2lmSm "C 1. 

An exact treatment in the framework of QFT with unstable vacuum (see, for example, [H [551 150] ) 
yields the following expressions for the scattering P{—\—)m of a particle (and an antiparticle) and pair- 
creation P(H — |0)„i probabilities, respectively (see Subsection 2.1 for notation): 

Pi-\-)m = \wi-HrnJP\ P(+ - |0),„ = \wi+-\0)^JP\ (66) 

where, due to relations ([T8|), (|19p and ([T6)) . the corresponding relative probabilities are 

k(-|"UJ' = YT^' k(+-|o)„„j' = Y^^- (67) 

As long as the semiclassical approximation is concerned (P" 1, H,„ <C 1), we have 

P(+-|0),„«h(+-|0)„„|'«H™. 

Thus, we can see that the quantities P(H — |0)m, \w (-1 — |0)mml^ ^'^^ ^"i '^^^^ identified only in the 
approximation under consideration. An exact expression for P" in terms of the mean values follows 
from ([5|), ([07)1 and reads 



P^ exp I K ^ In (1 - kH„) I 



(68) 

Formulas (|60p and ([05)1 imply the following exact relation between ImSm and H^: 

2ImS'™ = -Kln(l - kH„) . (69) 
It has to be used in the general case when the WKB approximation is not applicable. 
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3.3 Is it really a thermal distribution? 

Considerable attention has been recently focused on a mechanism of fast thermalization in heavy-ion 
collisions (see [TH |31] and references therein) . A possibility is discussed of a thermal one-particle dis- 
tribution due to quantum creation of particles from vacuum by strong electric-like fields. Some of these 
distributions are known in QED, and their relation to thermal spectrum of Hawking's radiation has been 
discussed (see references in the next section) . We give some examples of such distributions in Subsection 
4.3. However, a thermal one-particle distribution of created particles does not guarantee the character 
of thermal equilibrium for the corresponding complete quantum state of the system and only mimics, in 
some sense, the latter state. 

One ought to say that, in contrast to the case of Hawking's radiation, in which we do not have 
any information about another member of each created pair behind the horizon, both particles and 
antiparticles created from vacuum by chromoelectric field can, in principle, be observed. Because of the 
one-to-one correspondence between one-particle correlation functions and the reduced density operator 
of the complete system, all the moments of particle (or antiparticle) distribution coincide. However, the 
higher moments of the simultaneous distributions of particles and antiparticles are different. 

In what follows, we present a formal analysis of the above problems. 

Suppose that the differential mean numbers Nm of particles/antiparticles at the final state of a system 
subject to an external field have the form of a one-particle thermal distribution. There arises the question 
if one can be sure if in such a case the final state of the complete system is in thermal equilibrium, or the 
thermal form of a one-particle distribution can appear even if the final state of the complete system is 
not in thermal equilibrium. To answer these questions, we plan to examine two different possibilities of 
having the same one-particle thermal distribution for two distinct states of the complete system, one of 
them being a thermal equilibrium and the other a pure state. Let the first state of the complete system 
be described by the thermal density operator: 

= 1 exp {- [atij(+)a + fot } , 
Z = exp JK^ln(l + Ke-^"^') In (l + Ke"^-"') I , (70) 

n m ) 

where E'^^^ are given by (j27p . It is obvious that in such a state the differential mean numbers Nm have 
the form 

iV„- (e^" (71) 

On the other hand, if we have a causal evolution from vacuum, the density operator of the corresponding 
pure state having the form ((24|) : see the normal form in p7|) . Such a state provides the differential 
mean numbers ([7T|) in case (|57| holds true. We can see that measuring the one-particle distribution 
cannot distinguish between both these cases. Nevertheless, they can be distinguished by measuring the 
next moments, as demonstrated below. Let us calculate the variances Var^ in the states described by 
the density matrices ([70)) and ([^ . respectively, 



Vari'^ = tr 



Var^^ = tr {al^am + blnbm - 2iV™)^ p„ 
Since the differential mean values coincide in both states, one can see that 

Var^^ - Var^, = 2 (Q*„'^ - Ql^) , 

Qm = tr [alarr^blbmPp"'] , Qm = tr [at„a,„6,J„6™p^] . (72) 

To calculate the quantities Q^, and demonstrate that they are actually different, we are going to 
use the generating functional (J) ([SO]) , and the generating functional of momenta for the thermal 
distribution, 

^^1^ (J) = tr {exp pj(+)a + fe^J^-^fe] pf'^ . 

Then 



= K , (73) 



j=o 
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where the expression 
is used. On the other hand 



^exp|K^hi(^ 



I + Ke 



-£(+) + ,/(+) 



:^ln(l 



Nm [1 + (1 - k) N„ 



(74) 



.7=0 



Thus, we can see that the variances of the simuhaneous distributions of particles and antiparticles, Var^ 
and VarJ^, are quite different: 

Var*; - Var;; = 2iV™ {kN^ ~ 1) • 



4 Particle-creation in an electric-like background 



4.1 Quasi-constant electric field 

Below, we consider a number of applications of the above-developed formalism in QED with a quasi- 
constant (slowly varying) uniform electric field violating the stability of vacuum. We emphasize that our 
consideration can be relevant in QCD with an electric-like colour field and in some QFT models with 
a curved space-time, as was demonstrated, for example, in I40j. It was shown [J that the distribution 
of pairs created from vacuum by a quasi-constant electric field has a thermal-like form. It appears that 
such a form has a universal character, i.e., it emerges in any theory with quasi-constant external fields; 
and when applied to particle-creation in external constant gravitational fields it reproduces exactly the 
Hawking temperature. Thus, our consideration of QED with a quasi-constant electric field allows one to 
reveal the typical properties of any strong-field QFT. 

Note that in the case under consideration particle states are specified by continuous quantum numbers 
of the momentum p and spin projections r = ±1 (we formally set r = for scalar particles). From now 
on, we suppose that the standard volume regularization is used, so that (5(p — p') is replaced by (Jp.p' 
in the normalization conditions. Thus, our particles are labeled by a set of discrete quantum numbers, 
m = (p,r). 

As usual, we shall describe the electric field by time-dependent vector potentials. The states of the 
quantum system in question are far from equilibrium due to the field influence. We shall study in detail 
the time dependence of various mean values, in particular, the mean values of created particles. In a 
physically correct setting of the problem, we consider a model of a quasi-constant electric field E(x^) 
which effectively acts only for a finite period of time T and is zero outside this interval (we further call 
it the T-constant field). In our model, E(x^) = E for h < x° < ^2, ^2 = ~ti = T/2. Thus, the field 
produces finite work in a finite space volume. We accept the initial vacuum to be a free-particle vacuum. 
A relevant calculation in QED with a T-constant field can be found in [4]. Below, we use these results 
for evaluating the leading terms in particle-creation phenomena at large T, when the effects of switching 
on and off are negligible. 

Let us describe the T-constant field. It is nonstationary but has a constant space direction. We place 
the latter field along the cc^-axis. We denote by q the charge of a particle (by —q that of an antiparticle), 
and by M we denote mass. The corresponding potentials can be chosen in the form: Aq — Ai — A2 — 0, 
and 

r Eh, 

A3{x')^{ Ex\ 

[ Et2, 

where the time intervals are / = (— oo,ti), // ~ [^1,^2] 
If the time T is sufficiently large. 



x" e /// . 

///= (i2,+oo). 



T»To^il + X)/V\qE\, 



the differential mean numbers N„, read 




-VW\i<i< 



(75) 
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where K is a. sufficiently large arbitrary constant, K >> 1 + A, 



M^+pj 1 ^ \p3\-\qE\T/2 

P. = (^^-P,0), ^ . (76) 

and p3 is a longitudinal momentum of a particle [4] . One can examine the limit T — > cx) at any given p 
in the above expression. In such a limit, the differential mean numbers have a simple form: 

- e-"^ (77) 

which coincides with that obtained in a constant electric field by Nikishov [3]. One can see that the 
stabilization of the differential mean numbers to the asymptotic form ()77|) for finite longitudinal momenta 
is reached at T >> Tq- The characteristic time Tq is called the stabilization time. 

In order to investigate the effects of switching on and off for T >> Tq, we are going to consider a 
different example of a quasi-constant electric field: 

S(x°) = Scosh-2 (^^^ . (78) 

This field switches on and off adiabatically as ztoo and is quasi-constant at finite times. It is called 

an adiabatic field. The differential mean numbers of particles created by such a field have been found in 
|41j . For a further discussion, we need these numbers for a large a. As has been demonstrated in [Ij, the 
differential mean numbers in the field ([75)1 take the asymptotic form (|77p for a >> = {1 + VX) / y/\qE\ 
and for jpsj << \qE\ a. Thus, ao can be interpreted as the stabilization time for an adiabatic field. At 
the same time, the latter fact means that the effects of switching on and off are not essential at large 
times and finite longitudinal momenta for both fields. Extrapolating this conclusion, one can suppose 
that particle-creation effects in any electric field, that is one being quasi-constant sa i? at least for a 
time period T >> Tq and switching on and off outside this period arbitrarily, do not depend on the 
details of switching on and off. Thus, our calculations in a T-constant field are typical for a large class 
of quasi-constant electric fields. 

It is of interest for phenomenological applications to calculate the distribution of particles created 
with all possible values and a given p± (it is called the p± distribution and is denoted by Up^^ , in what 
follows). Analysing the total mean number of particles created by the T-constant field, we go over from 
summation to integration: (^Tp" / "^P' Then, the total mean number (we denote it by H) can be 

presented as 

\i = V J (fp^np^, (79) 

where 



1 



(27r)= 



jdP3^rn (80) 



is the p± distribution density of particles created per unit volume. Hm is constant for jpsj < \/\qE\ I y/\qE\T/2 — K 



and for T » Tq, and decreases rapidly for \p^\ > y'\qE\{^J\qE\T/2 + K\. The contribution to the 



integral ([80]) from the intermediate region can be estimated as 2^J\qE\K . This implies 



./\^\Te-^^ + 0{K) , (81) 



where J is the number of the spin degrees of freedom ( J = 1 for scalar particles and J ~ 2 for fermions). 
Thus, the pj_ distribution density of the particle production rate has the form 



dup^ J I qE I 



dT (27r)3 

The total number of particles created per unit volume is given by 



e"'^^ . (82) 



H (qE? T r M2 1 

V-'-^-A-'W\]^ (83) 
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Suppose that, in addition to the external electric field, there exists a parallel constant magnetic 
field B. For definiteness, let us choose its potentials as Aj^ = Bx^S^^. Then, the complete set of 
quantum numbers that describes particles in such a background is (pi, 7iB,p3, r), ub = 0,1,..., and 
A = (M^ + \qB\{2nB + 1 - r)) \qE\~'^ . Substituting this A into eqs. jTZ]) and we obtain the differ- 
ential mean numbers and p^-distribution density of particles created per unit volume, respectively. One 
can see that the presence of the magnetic field essentially changes the energy spectrum of the transversal 
modes (which used to be p^^ in the absence of the magnetic field) and also the form of A, but does not 
change the dependence of and Up^aiv A. This is related to the fact that the magnetic field itself does 
not produce work acting on charged particles, and, therefore, does not create particles. Then, the total 
number of particles created per unit volume is given by 

H Jq^EBT\cosh{TTB/E)]^^+'^^/'^ { AP ] , , 

- L V / ;j .exp -TT— L (84) 



V ^T:'^smh{TT B/E) \ \qE\\ 

where k = +1 for fcrmions and k = — 1 for scalar particles. 

For a strong electric field, IvP / \qE\ < 1, (B = 0), and large T, the energy density of created pairs 
reads £ — \qE\ TH/y; see |42]. We can neglect the back- reaction of particles created by the electric field 
in case their energy density is essentially smaller than the energy density of the electric field, £ <C S^/Stt. 
Consequently, the concept of a strong constant electric field is consistent only if the following condition 
holds true: 

l«ki.|T^«-^exp|.— |. (85) 
Following [4], we represent the asymptotic formula (j77p in a universal form: 

H„ = exp(-27r^l , (86) 



9 

where is the work of an external field creating a particle from a pair in a given state m, 

where Po{tf) and Po{ti) are particle energies at the final time instant tf and at an initial time instant 
ti, respectively, and Ae^ac is a shift of the vacuum energy due to the time evolution. The quantity g is 
the classical acceleration of a particle at the final time instant. In the case of a T-constant field, one can 
obtain 



M^ + pI _ X _ \qE\ _ 2 



2po{tf) T' ^ poitf) T 

Thus, we can see that the differential mean values (f86| are given, in fact, by the Boltzmann formula with 

a 



the temperature 6 = „ ^ (where ks is the Boltzmann constant), the latter having literally the Hawking 



form [8]; see below. 

We recall that the Hawking result for bosons created by the static gravitational field of a black hole 
in a specific thermal environment has the Planck form 



exp <2tt- 



9{H) 



(87) 



Here, ujm is the energy of a created particle and the Hawking temperature reads 9{^h) = where 
g{H) = is the free-fall acceleration at the gravitational radius rg of a black hole with mass M. In 
this case of a quasi-static gravitation field, the evolution shift of the vacuum energy is l^Cvac — 0, so that 
one identifies the work ojm (we have introduced) with the energy of a particle in formula (|87|) . It is also 
known [43^ that an observer which moves with a constant acceleration g^^) (with respect to its proper 
time) will probably register in the Minkowski vacuum some particles (Rindler particles) . The distribution 
of Rindler bosons has the same Planck form (|87p , where one needws to replace g(H) by 9{b) : so that the 
corresponding temperature is 6'(ij) = ^.tiVb ' 

It is a direct consequence of the equivalence principle that the effective temperature 9 of distribution 
has literally the Hawking form. The different form of distributions can be caused by essentially 
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different structures of the Fock space in both cases. We beheve that the Planck distribution arises neces- 
sarily due to the appearance of an event horizon (there is a boundary of the domain of the Hamiltonian) , 
that is, due to the condition for which the space domains of the particle and antiparticle vacua are not the 
same. On the other hand, the final state can be treated as an equilibrium state. In contrast to this, in a 
uniform electric field we deal, in fact, with both the particle vacuum and the antiparticle vacuum defined 
over the entire space, that is, these space domains coincide. In this case, the mixed state of particles 
(antiparticles) described by the iPv-) density matrix of Sec. 3 can be represented as a pure state 
in an extended phase space where the space domains of both the particle vacuum and the antiparticle 
vacuum are the same, being a state of a far-from-equilibrium system. Let us note that in the framework 
of a semi-classical description at oJmld << 1 the Boltzmann spectrum closely approaches the Planck 
spectrum. 



4.2 Soft parton production by SU(3) chromoelectric field 

As mentioned in Introduction, in QCD there exist physical situations that are quite efficiently described 
by the chromoelectric flux tube model. In this model, the back-reaction of created pairs induces a gluon 
mean field and plasma oscillations (see |^ and references therein). It appears that the calculation 
of particle-creation in this model requires the application of the general formalism of QFT for pair- 
production at a finite temperature and at zero temperature both from vacuum and from many-particle 
states (see [Ml ESI |45] for physical reasons) . The study of various time scales in heavy-ion collisions shows 
that the stabilization time Tq is far smaller than the period of plasma and mean-field oscillations. Then, 
according to condition (j85p , the approximation of a strong T-constant chromoelectric field can be used in 
treating such collisions during a period when the produced partons can be considered as weakly coupled 
due to the property of asymptotic freedom in QCD. It may also be reasonable to neglect dynamical back- 
reaction effects and to consider only pair-production from vacuum by a constant SU (3) chromoelectric 
field. 

Here, we would like to turn our attention to results obtained in QCD with a constant S'[/(3) chro- 
moelectric field E"" [a = 1,...,8) along the z'^-axis; see [12]. In this work, the imaginary parts of 
one-loop effective actions for quarks 5''"'"''^ and gluons 5"^'"°" have been calculated via gauge-invariant 

respectively. They have the form 



Pj^ distributions 



and 59^"°" 



uark 



ImS"^ 



ImS'^'^""" 



j (fpi_ Im S'«';'"^'= , Im 59'"°" = J (fp^ Im , 



VT 



VT 
8^ 



In 1 + e 



i)l 



Hi) 



Pi 



where £'(j') are the eigenvalues of the matrix iT'^E'^ for the fundamental representation of SU{3); E(^j-^ 
are the positive eigenvalues of the matrix if°-^^E^ for the adjoint representation of SU{?>); and q is the 
coupling constant. These eigenvalues are the following gauge-invariant quantities: 



^(1) === v/C'i/3 cos 61, E(^2) ^ v/C'i/3cos(27r/3- 
E(^) = v/c73cos(27r/3 + 0) , 

where 9 is given by cos^ W — 36*2 /Cf, and 



Ell 



(1) 



E, 



(3) 



— 1 — cos 
2 V 

— 1 + cos , 
2 V V3 



1/2 



' ^(2) = 
1/2 



1 + COS 



f- 

V3 



-,1/2 
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where 6 is given by cos^ ^ = — 1 + 6C2/C1. Here, Ci and C2 are Casimir invariants for SU{3), 

where dabc is a symmetric invariant tensor in the adjoint representation of SU{3). Then, the probabihties 
for a vacuum to remain a vacuum are found, for both quarks and gluons, from relation (j60p . However, 
formulas for parton production rates obtained in [12] hold only in the approximation 2 Im S'^""'''' <C 1 and 
2ImS'|'^"°" ^ 1 by virtue of the arguments that we present at the end of Sbsection 3.2. To obtain exact 
results, we can use the following line of reasoning. The results ([88]) can be treated as those obtained in 
the case of a T-constant chromoelectric field when the integration over the longitudinal momentum and 
the summation over the spin and color degrees of freedom have been carried out. Then, using relation 
(|69p , we can extract from representation (|88p an exact expression for p± distribution densities of quarks 
^guarfc ^^^^ gluons n^*^"" produccd per unit volume. Those are 



e^'^^o) , (89) 



where T is a sufficiently large action period of a constant field. The pj_ distribution densities of particle 
production rates can be found as dn'^^^'' / dT and dn^''^°^ /dT, respectively. The total numbers of quarks 
and gluons created per unit volume can be obtained from ([55]) as follows: 

^quark 7^ 3 r M'^ ] 3Tq^Ci 

where the relation J2^=i E'^j) — 3Ci/2 is used. Taking into account the relation -^fj) ^ one 

can see from ([90| that in a sufficiently strong field E(^j) , / | ^ 1, the densities of created quarks 

and gluons are related by '^i^'^^k jy ^ •^giuoni^y^ 

We can see that the {])- terms in expressions (|88)) and (|90p can be interpreted as those which are 
obtained for Abelian-like electric fields -Eq) and , respectively. The maxima of the fields are restricted 

by the conditions < Ci/3 and Eq) < ^fCi. Therefore, in order to study the validity of the 

constant SU{i) chromoelectric field approximation we need to take into account only the energy density 
of gluons created by the field -Eq). We know from the previous subsection that this energy per a single 



pair IS 



T. Then, the total energy density of created gluons reads 

< lol JCiT- 



3 3 ^gi'^ 



T 

T^E 9^0- 



47r3 ^ 



V 



One can neglect the back-reaction of these gluons created by the chromoelectric field only if £ <cCi/87r. 
Finally, the condition of validity of the T-constant SU{'i) chromoelectric field approximation can be 
written as ^ 

1 « kl \fclT^ « ^ . (91) 



Therefore, we can see that the T-constant SU{i) chromoelectric field approximation is consistent 
during the period when the produced partons can be treated as weakly coupled. 

Recently, it has been discovered (see, [46]), in the framework of the Color Glass Condensate Approach 
to the description of heavy ion collisions, that shortly after the collision the system contains both longitu- 
dinal chromoelectric and chromomagnetic fields. In such a chromomagnctic field, the above-described soft 
parton production is changed. First of all, the energy spectrum of the transversal modes becomes differ- 
ent, which implies that the expressions for that enter Aq) and A(j) are different from the case of pure 
chromoelectric field. For example, if the chromomagnetic field B'^ is uniform then = \qB'^\{2nB + ^~r), 
ub = 0, 1, . . .. This expression is similar to the case of QED, see Section 4.1. Here, the spin quantum 
number r takes the values ±1 for quarks and —2, 0, +2 for gluons. As in QED, the constant chromomag- 
netic field does not produce work acting on charged particles, and, therefore, does not create particles. 
The distribution densities of created quarks n^""'''^ and gluons functions of Aq) and Aq) do not 

change and are given by equations 
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Nevertheless, the above expression for in a uniform chromomagnetic field is negative for ub = 0, 
r — +2. This is why such a configuration of the chromomagnetic field is unstable under quantum 
fluctuations. For the first time, it was mentioned in [47| . Then, the problem has been discussed in 
numerous articles: for a review, see, e.g., [48l H^. One of the possible stable configurations of the 
longitudinal chromomagnetic field is the so-called flux-tube state, alias the spaghetti state (see [l^ and 
references therein). An explicit form of the spectra p\ in this field is unknown. Nevertheless, a stable 
chromomagnetic field is constant and longitudinal. This is sufficient for us to make the above conclusion 
that the dependence on A(j) and Aq) , given by eq. (j89p . in the presence of a longitudinal chromomagnetic 
field remains unchanged. 

In the following sections, we turn once again to particle creation by electric field in QED. The above 
discussion shows that it can be useful for understanding the effects of quark and gluon creation in QCD. 



4.3 Thermal-like distributions 



As has been mentioned, the thermalization stage of multiparticle production in ion-ion collisions at high 
energies is very important. On the other hand, as we know from Section 3, it is sometimes difficult to 
distinguish a real thermal equilibrium from a state where we have a one-particle thermal distribution. In 
this connection, we shall consider some simple examples when pair-creation by electric field can mimic a 
one-particle thermal distribution. 

We recall that due to the screen of created pairs, the original electric field may have an exponential 
fall-off: 

E{x^) = £;e-^"/" . (92) 

The differential mean number of particles created from vacuum by this field has been calculated in [50| . 
The result is 



cosh 


7rQ(£4-P3) 
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™(e-P3) 


■e 


sinh 


7rQ(e+P3) 




sinh 


7I-q(£-P0 


-e' 



1 



Bose case, 



Fermi case, 



(93) 



where e = y^M^ + j? and — p3Sgn{qE) . For na jpaj ^ 1, these expressions coincide with the Bose 
and Fermi distributions at the temperature 9 = {2TTkBa) ^, respectively. 

Another example is the pulse of electric field (|78p . The differential mean number of particles created 
from vacuum by the sharp field pulse ([75]) at a \ qE\ /e <C 1 can be extracted from the result [31] and has 
the form 



{nqEa^y 



{qEa^)U 
1 - {Pz/e) 



(ps/s) sinh {nae) Bose case, 
sinh^^ (nae) Fermi case, 



(94) 



When the ratio |p3 1 /e is sufficiently small and the external field is not strong, e ^ \/\qE\ ^ 1 , there is a 
range of values a, nae 3> 1, in which distributions (j94p have the Boltzmann form with the temperature 
9 = {2'KkBa)~\ 

As we note in the previous subsection, depending on the details of the chromoelectric fiux tube model, 
its stage and field strength, pair production both from vacuum and from many-particle states by a T- 
constant electric field may be relevant. The well-known asymptotic form (|77p of the differential mean 
numbers for pairs created from vacuum by a constant electric field can lead to a thermal-like distribution 
of created pairs if, for example, the chromoelectric string tension undergoes Gaussian fluctuations |51j . 
This implies a modification of the original flux tube model by introducing a fluctuating string tension. 
In case the original flux tube model still holds, nevertheless, the differential mean numbers of pairs 
created by a T-constant electric field can be represented as the Boltzmann distribution (|86p with the 
temperature 9 — (jrkBT) , as has been seen in Subsection 4.1. Then, it may be reasonable to examine 
the phenomenological model with a slowly oscillated mean electric field and suppose that the pair creation 
by the mean field during the semiperiod of oscillation can be effectively approximated by a T-constant 
electric field (we recall that the time scale of stabilization Tq is far smaller than the period of oscillations). 
In this case, the electric field produces pairs for a semi-period of oscillation in the presence of pairs created 
at previous stages. This is the way to take into account the effects of back-reaction in such a model. 
In other words, we are going to consider pair-creation from the initial state given by a distribution of 
previously created particles. Formula (|59|) is relevant in this analysis. 
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Starting from the initial vacuum state, one has AA^^ = ^m, where belongs to the asymptotic 
form ([77]) . Then, at the end of the first stage, when the mean field is depleted for the first time, the 
distribution of particles (being equal to that of antiparticles) is Nm'^ — N„j . During the second stage, the 
direction of the mean field is opposite to the field direction at the first stage. Due to the condition of 
stabilization, this is of no importance, since is an even function of qE. Thus, when the mean field is 
depleted for the second time, from ([SS]) there follows the equality 



and at the end of the n-th stage, 



Consequently, the total number of particles created at the end of the n-th stage is 



n-l 



N^a^ = V (1 - 2«;N„)' 



i=0 

We have this result if the created particles do not leave the region of the active field. In order to take 
into account the possible loss of particles due to interaction, movement, etc., we also assume that the 
total number of particles at the initial state of the n-th stage is less than the number N^^~'^^ of particles 
created at the end of the (n — l)-th stage and is 7A^^"~^\ where 7 < 1 is the factor of loss. Then, the 
modified relation is 

TVi,"' = H„ (1 - 2k^^) 7A^i"-'' , (95) 

and we finally have 

n-l 

^^r^ = N„ ^ 7' (1 - 2K^^r,^)' ■ (96) 

Supposing that 7 is constant, one can calculate the sum in (|96p . 

^i"^ = ^m- , r = 7 (1 - 2^cH™) . (97) 

1 — r 

For fermions, k, = +1, then Nm^ < 1. Energy dissipation after a period of oscillation is estimated 
(for real parameters of heavy-ion collisions) not to be large, so that damping is small and the number 
of oscillations can be quite large; damping decreases with an increasing field strength. If the number of 
cycles is sufficiently large, we get the limiting thermal-like distribution 

^ ^ ^ i I f98) 

™ 1-7(1-2kN„0 27'e'^^(l-7)/27 + K' ^ ' 

In other words, the system reaches a quasi-equilibrium state. For bosons, k = — 1, then Nm"^ increases. 
This is the phenomenon of resonance, and the increase can be either limited or unlimited depending on 
the factor 7. The increase is limited as long as r < 1. In this case, formula (|98p is valid for bosons as well. 
We can see that back-reaction-induced plasma oscillations can reach a quasi-stationary form specified by 
thermal-like distribution for both bosons and fermions. 



4.4 Particle creation at finite temperature 

We are now ready to present explicitly the mean number of (anti) particles in the mode m (with finite 
longitudinal momenta, jpaj < \/\qE\ (^^/\qE\T/2 — ) for the final state of evolution in a quasi-constant 

field from the initial thcrmodynamical equilibrium, Nm\in) = (e^™ + k) ^, at equal chemical potentials 
= /i'^' ~ {11 < M for bosons), 

Ar(C) ^ (e^^™ e--^ (tanh {E,r./2)f , (99) 

where Em = P {£m — m); = ^ M'^ + pi + (tts)^, t^s = Ps + qET/2, and it is implied that A is given 
by l|75|) . This result for the electric field coincides with the one obtained in [M]. Due to the effect of 
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stabilization, it seems that the time-dependence of the final distributions in question is absent. However, 
the integral mean numbers vary as long as a quasi-constant field is active. 

It is of interest to establish a general behavior of the integral mean numbers of created particles when 
the effects of switching on and off are negligible. As has been shown above, we can satisfy this condition 
by selecting the action time T of the T-constant field (T >> Tq) as an effective period of pair creation. 
It is implied that, in general, the final time instant, tf, and the initial time instant, ti, are so selected 
that the quasi-constant field is closely approximated by the T-constant field for a period from ti to tf, 
and tf — U = T. 

Let us estimate the sum over the longitudinal momentum of ANm in (|59|) , which is the mean number 

■ / dp, and the distribution 



of particles created with all the possible values p^. As above, 

plays the role of the cut-off factor for the integral over p^. Then, one can conclude that the p±,r 
distribution density of particles produced per unit volume is finite and can be presented as follows: 



1 



{2n) 



AN„,dp3 



\qE\T/2 



n^iP)dp3 + x/\^\OiK) 



-\qE\T/2 



n™(/3) = (tanh(i;™/2))'' 



(100) 



From (|100p , one can estimate the p± 

dn'i 



r distribution density of the particle production rate, 
19^1 ^ ..M 



dT 



"m(/5)U3 = | 



qE\T ' 



(101) 



We have [qET) ^ M + p]^, according to the condition of stabilization. Then, the high and low 
temperature limits for the production rate are only defined by the final longitudinal kinetic momentum 
\qE\ T and the temperature 8 relation, /? \qE\ T ^ 1 and /3 \ qE\ T ^ 1, respectively. For simplicity, we 
assume that \qE\ T ^ jj,. Considering these limits, one obtains for the temperature-dependent term in 
(|T0T|l 

^^nm.,=\qE\T = l-2Ke-'3|*^l^, /3|g£;|r»l, 
nmm.,=iqE\T = [P\qE\T/2r, P\qE\T^l. 



We can see that at high temperatures the rate 



dT 



is time-dependent; it is much lower than the zero- 



temperature value but increases for fermions, and is considerably higher than the zero temperature value 
but decreases for bosons. Consequently, the frequently used notion of a number of particles created per 
unit of time makes sense only at low temperatures and in this limit it coincides with the zero-temperature 
value of the production rate. We consider two temperature limits for the p±, r distribution density (jlOOp : 
low temperatures at f3 {e± — /i) ^ 1, = ^/M^ + p^, when all the energies of the particles created in 
the modes with a given p^ are considerably higher than the temperature O, and high temperatures, at 
f3\qE\T <^ 1, when all the energies of the created particles are much lower than the temperature 8, 



P±,r 



{2^f 
(i\qE\ 
{2nf 



Pi,!" 



VM^Te-'^^ + O (if ) ,K-±1, /3(e_L-M)>l, 



{2^f 



\qE\T^/2 + 0[.J\^\T 
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-1, /?|g£;|T< 1, 



(102) 



\n[./\^\T/K)e-''^ + 0{K) 



, K= -1, l3\qE\T<^l. 



The result at low temperatures is not different from the zero-temperature result [1] within the accuracy 
of our analysis. Integrating expressions (|102[) over pj^, one finds the total number of particles created per 
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unit volume at the low-temperature and high-temperature limits, respectively: 



TV 



V " (2^)3 



^cr \qE\\n{^^/\^\ 



T 



V 27r3/3 



e-"^^'/^^, K = /3|g£;|r< 1, (103) 



where the summation over r = ±1 is carried out for the fermions, and only the leading T-dependent terms 
are presented. From (|102p . (|103p . we can see that the values of the integral mean numbers for fermions 
at high temperatures are much lower than the corresponding values at low temperatures. For bosons, 
the integral mean numbers at high temperatures are considerably higher than the corresponding values 
at low temperatures. 

As mentioned in Introduction, thermally-influenced pair production by a constant electric field has 
been investigated in several approaches PB] [551 IH] ■ The results are quite contradictory, varying from 
the absence of creation to the rates of fermion production higher than the rate at zero temperature. Now, 
we are ready to discuss these contradictions. As has been shown above, the initial thermal distribution 
affects the number of states in which pairs are created by a quasi-constant field. Hence, pair-production 
exists at any temperatures, and, in particular, the fermion production rate cannot be higher than the 
rate at zero temperature, by any means. Note that our calculations are based on the generalized Furry 
representation developed especially for the case of vacuum instability in accordance with the basic prin- 
ciples of quantum field theory. On the other hand, all the conclusions of 26, 27 concerning the pair 
production rate and/or the mean numbers of pairs created at non-zero temperatures are based on either 
the standard real-time, or imaginary-time, one-loop effective actions. However, such formalisms do not 
work in the presence of unstable modes. The real part of the standard effective action describes the 
effects of vacuum polarization and has nothing to do with the time-dependent conduction current of 
created pairs. For example, this can be observed at zero temperature (see [42]). In this case, the in- 
formation concerning pair-creation comes from the imaginary part of the standard effective action. The 
extension of real-time techniques for finite-temperature quantum electrodynamics with unstable vacuum 
has been presented in [32] ■ In this article, one can see that the relevant Green functions in a constant 
electric field are quite different from the standard proper-time representation given by Schwinger. Then, 
the relevant real-time one-loop effective action must be different from the standard one[£|. The standard 
imaginary-time formalism, obtained under the assumption of a thermal equilibrium and the appearance 
of a contradiction with the Pauli exclusion principle, shows that the attempts of generalization to far- 
from-equilibrium systems have failed. The functional Schrodinger picture used in [25] to calculate the 
TV°'" at high temperatures seems relevant; its asymptotic expressions for TV'^'" at high temperatures are in 
agreement with our expressions in (|103p . 
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Appendix 

I. For both the Bose and Fermi cases, the following relations hold: 

^gatDa ^ ga+DagD^ ^ a) e''^ = e"'' ^ " , (104) 

e'^'^'^=:exp{at(e^-l)a} :, (105) 



^We will present the relevant real-time one-loop effective action elsewhere. 
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where D is a matrix. To prove (jlOSp . let us consider the operator function F{s) = e*°*-^°, where s is a 
parameter. The function is a solution of the following equation: 

^^^a^DaF(s), F(0) - 1. 
as 

Using relation p04p . we can rewrite the right-hand side of the equation as follows: 

^^^a^F{s)De'°a, F(0) = 1 . 
as 

Now, we can verify that a solution of such an equation reads 

F(s)=:exp{at(e«^-l)a} : . 

Setting s = 1, we justify (|105p . 

II. We often use the well-known relation 

where A and A arc Grassmann-odd or Grassmann-even variables depending on statistics. For a product 
of two normal forms, there holds a generalization of (|106p . namely, 

• 6°^-°" • ■ e"^^" ■ = ■ e°-^(D+D+DD)a . (107) 

where D and D are some matrices. 

III. The projection operator on the vacuum state can be written as follows: 

Po = |0>(0| =: e-''^" : . (108) 
Such a representation was first used by Berezin [33J . One can see that the operator Pq obeys the equations 

aPo =0, Poa^ -0, Po|0) - |0) . 

Using the Wick theorem, one can see that : e""^^" : is a solution of these equations. 

IV. The trace of a normal product of creation and annihilation operators can be calculated by using 
the following path integral representation. Let X (a^ , a) be an operator expression of creation and 
annihilation operators, a and a^. Then the trace of its normal form 

oo 

tr{: X{a\a) :} - ^ ^ (A/!)-^0|a™,, . . . a™, : X (at, a) : a^, • ■ • a^^jO) , 

M=0 {m} 

can be expressed as the following vacuum mean value: 

tr{:X(at,a) :} = {0\T : X {a^ , a) : e''(*^)'^'(*')|0) (109) 

where the notation a = a {tf), = (ti) is used for the operators a to the left of : X (a^, a) : and to 
the right of : X (a^, a) :, whereas T is the ordering operator putting a{tf) to the left of : X (at, a) : and 
at (ti) to the right of : X (at, a) :. Using either the Berezin path integral or the Gaussian integral over 
ordinary variables, depending on statistics, one can rewrite (|109p as follows: 

tr{:X(at,a) :} = (0| J exp {kA*A + A*a} : X (at, a) : exp {atA} ndA*dA|0), (110) 

where a (tf) — a and at (t,) = at are used after rewriting. 
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